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EQUATION OF MOTION FOR DISTRIBUTED SYSTEMS 


Differential equation: Lu(P^t) + M(P)3^u(P^t)/Jt^ = f(P>t) ^ Pe D 

Boundary conditions: Bju(P>t) = 0^ l=L2>...^p ^ Pes 

u(Pyt) = displacement at point P 

Ij Bj = differential operators (L is self-adjoint of order 2p) 

M = mass density 

f(P^t) = distributed control force 


EIGENVALUE PROBLEM 

Differential equation: L4*=AM+ 

Boundary conditions: Bj«#»= 0, i = L2^...,p 
Solution: eigenvalues A ^ =o^ , eigenfunctions (r=L2.,,.) 
Because L is self-adjoint^ eigenfunctions are orthogonal 
L is generally positive semldefinite -► Ap are all nonnegative 

Ap = 0 for rigid-body modes 

= natural frequencies 

Orthonormality conditions: M<f»s+r clD = Jpg . fgLfp dD = Apfp^ 
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MODAL EQUATIONS 


Expansion theorem: u(P>t) = Z 

r=i I I 

Up(t) = modal coordinates 

Modal equations: Up(t) +o^Up(t) = fp(t) > r=L2>... 

Modal controls: fp(t) +p(P)f (Pyt)dD ^ r=h2^... 

Coupled controls: fp(t) = fpiu^^ Ui^ U2-' ■' f=1^2^... 

Independent modal-space control (IMSC): fp(t) = fp(Up^Up)^ r=L2^... 


CONTROL IMPLEMENTATION 

Distributed actuators: f(Pyt) = Z 

Coupled controls: unable to design distributed controls 

IMSC: design modal controls firsts then use above formula 

no control spillover 

m . 

Discrete actuators: f(P>t) = Z F|(t) AIP-P.) 

j = l J J 

m m 

f.(t) = Z = 21 BriFj(t) > r=L2..... 

r j=i T'r j j j=i rj j 

n = number of controlled modes 


i-[Pl F2 




f = 0, 


fjl B - B 


ri 


Coupled controls: Design f(t) so as to ensure controllability 
IMSC: Design f (t) so^as to control a given number of modes 

Then^ FTt) = B^f(t) . B^ = pseudo-inverse of B 
To avoid pseudo-inverses^ let m = n^ or the number of actuators 
must equal the number of controlled modes. 


CONTROL IMPLEMENTATION (CONT^D) 


Distributed measurements: 

Measurements: u(P^t)^ u(P/t) for all P and at any t 

Then^ modal coordinates and velocities^ Up(t) and Up(t)^ are computed 

by using the modol filters 

Ur(t) = ( H(P) +r(P)U(P.t)dD. Ur(t) = ( M(P) r=L2. . . . 

Discrete measurements: 

Measurements: yj(t) = u(Pj^t)^ yj(t) = u(Pj^t)^ j = h2>...A 
k = number of sensors 

Standard approach: Use Luenberger observer to estimate state 

Discrete measurements treated as distributed: 

Use interpolation functions to compute estimate ti(P/t) of u(P^t) 

Then^ use modal filters to compute estimates Uj-(t) of Up(t) 

Divide structure into s segments (elements) 

A ^ T 

Approximate displacement: u(P^t) = E L (P)Vi(t) 

j=l 

Vj = measurements at the boundaries of j'th interval 

Lj = vector of interpolation functions (from the finite element method) 

Estimated modal coordinates: 

s s 

Ur(t) = S M(P)+r(P) H L|vi(t)dD= E IriVi(t). 
r r -j~j -rj-j 

r = ly2y...yn 

c s 

Similarly Up(t) = E IriVj(t) 
r j = j ^ f J ^ J 

irj " lj^M(P)+p(P)Lj(P)dD = const 

are computed off-1 ine^ in advance. 
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CONTROL IMPLEMENTATION (CQNT^D) 


Rearrange Iri such that 

Up(t) = i Cfjyj(t) . Uf-(t) = ^rJ 

Let 

u(t) = [uj U2‘" ^2 • • ■ Up]^ 

C “ C p j J r~ 1 ^ 2 ^ ■ • » j rij j~ 1 ^ 2 ^ ■ • i^k 

y(t) = [yi y2 ... y(t) = [vi 1/2 .•■ yj^ 

u(t) = Cy(t) y u(t) = Cy(t) 

M 

The way C is assembled depends on the nature of the interpolation 
functions (see example later). 


THE LANGLEY BEAM EXPERIMENT 


Free-free beam controlled by using 4 actuators and 9 sensors 

Use IMSC to control four modes^ two rigid-body and two elastic modes 

Actuator forces Fj(t) = (B"^)f-jff(t)^ j = L2^3>4 

Modal forces: 

1) For rigid body modes 

7 p = I Upl +1 Upl/Cp . r=h2 
Cp = weighting factor 
if 7 p < dpy then fp = 0 
if 7 p > dp and 

i) Up > 0> Up > 0> or Up > 0 > Up and lUpl< fp> then fp = -kp 

i i ) Up < 0 ^ Up < 0 ^ or Up < 0 < Up and I Up)< c^j then f p = kp 
dp = magnitude of the deadband region 

£p = threshold velocity^ kp = modal control force 



THE LANGLEY BEAM EXPERIMENT (CONT'D) 


2) for elastic modes 

fj-(t) = -kj. > |Ufl » dp and Up > 0 

0 y -dp < Up < dp 

kp ^ / Up I ^ dp and Up < 0 

For simulation purposes^ the response is available in closed form both 
for rigid-body and elastic modes. 

Actuator locations immaterial when IMSC is used. 


F(t) = B"^f(t) Fj(t) = z: (B‘^irfr(t) j=L2,3.4 
~ ~ J r=l ' 


Actual control forces are a combination of on-off functions of the type 


-Jr 


fr 

Ur 

Jr 7r 

•Ur 


Actuators available at LARC have four components/ which can be assigned 
to each modal control force. As a result/ each actuator command becomes 
a linear combination of 4 modal on-off control forces. 



THE LANGLEY BEAN EXPERIMENT (CONT'D) 


Sensors measure displacements alone. Ideal ly> for bending they should 
measure displacements and slopes (velocities and slopes of velocities too) 

Divide the beam into four equal segments (elements) and measure displacements 
at the ends and middle point of the elements so that the nine sensors are spaced 
at equal intervals. 

As interpolation functions use 

Lj = [i(2^D l-3t+ 0< t< \ 

where ^ is a local coordinate related to the global coordinate x by 

i = j - in which j is the segment number 


The C matrix is assembled fromJ^pj tensor as 

^rp " Sl-rJi ^2j+l-2>p ^ ! U2,...jn , n=A 

p = \j2j k = 9 


£ = index denoting the interpolation function 
j = index denoting the element number 

Because velocity measurements are not available estimate velocities by using the 
relation 

Ur(jT) - Ur(JT-T) 

Up(jT) = 

T = sampling time = 1/33 sec. 

Or^ one could use a modal Luenberger type observer. Because the controls 
are nonlinear^ the convergence of the observer can only be determined by trial 
and error. 
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THE LANGLEY BEAM EXPERIMENT (CONT*D) 


Parameters associated with the beam: 

L = 12 ft. cross-section = 6 x 3/16 in 

6061 aluminum: p= 0.1 Ib/in^ E = 1 x 10^ Ib/in^ 


The free-free> 
transcendental 

Oj = 0 
^2 = 0 

0 J 3 = 11.47979 
= 31.64450 
W 5 = 62.03586 
<Jg = 102.5484 
«7 = 153.1897 


uniform beam admits a closed-form solution, 
equation was solved numerically to yield the 


rad/s 

rad/s 

rad/s 

rad/s 

rad/s 


The 

eigenvalues 


Simulation of the Beam motion : 

The first 7 modes are Included in the simulation: 4 controlled + 3 
residual modes 

Control Gain Parameters: 

dj = d 2 = 0.002 . d 3 = di( = 0.0005 

kj = l <2 = 0.3 . k 3 = 0.12. ki, = 0.03 

= £-2 = O'Ol 

= lUpl + lu^l /lO . r=1.2 

Sampling time = 1/33 sec. 

Viscous damping was added to each flexible mode 
Damping factor = 0.002 . r=3.4.5.6.7 

Disturbance of the beam was taken in the form of a unit Impulse of magnitude 1/12 lb 
applied at Xg = 0.67L. 

The displacement of the beam cannot exceed 1 in at any time because of the location of 
the sensors and actuators in the experimental setup. 


118 



. IN 















THE LANGLEY BEAM EXPERIMENT (CONT*D) 


Results 


• The main contribution to the response is from the rigid-body modes. 

•The second elastic mode shows noticeable participation. This participation will 
eventually disappear due to internal damping. 

• Control of the second elastic mode can be enhanced by sensing velocities^ or 
estimating velocitites via a Luenberger observer. 


• Observation spillover (which may arise from the need of more sensors) was found 
to be negligible. So was the control spillover into the residual modes. 
Simulations of the beam with and without the residual modes indicated that 
spillover effects are infinitesimal. The reasons for this are: 

1) IMSC is used 

2) Modal filters are used 

3) Residual modes have very high frequencies 


Conclusion 

The IMSC method in conjunction with on-off modal controls does a good Job 
in controlling the motion of the beam^ where the motion consists primarily of the 
rigid-body modes. 
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